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ABSTRACT 

We present our temporal and spectral analyses of 29 bursts from SGR J0501+4516, detected 
with the Gamma-ray Burst Monitor onboard the Fermi Gamma-ray Space Telescope during the 
13 days of the source activation in 2008 (August 22 to September 3). We find that the Tgo 
durations of the bursts can be fit with a log-normal distribution with a mean value of ^ 123 ms. 
We also estimate for the first time event durations of Soft Gamma Repeater (SGR) bursts in 
photon space (i.e., using their deconvolved spectra) and find that these are very similar to the TqqS 
estimated in count space (following a log- normal distribution with a mean value of ~ 124 ms). We 
fit the time-integrated spectra for each burst and the time-resolved spectra of the five brightest 
bursts with several models. We find that a single power law with an exponential cutoff model 
fits all 29 bursts well, while 18 of the events can also be fit with two black body functions. We 
expand on the physical interpretation of these two models and we compare their parameters and 
discuss their evolution. We show that the time-integrated and time-resolved spectra reveal that 
-Epcak decreases with energy flux (and fluence) to a minimum of ^ 30keV at F = 8.7 x 10^^ 
erg cm~^ s~^, increasing steadily afterwards. Two more sources exhibit a similar trend: SGRs 
J1550 — 5418 and 1806 — 20. The isotropic luminosity, Liso, corresponding to these fiux values is 
roughly similar for all sources (0.4 — 1.5 x 10'"' erg s~^). 

Subject headings: soft gamma repeater: general — soft gamma repeater: individual(SGR J0501-f4516) 
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1. Introduction 

Magnetars are slowly rotating neutron stars as- 
sociated with extreme magnetic fields {B > 10^^ 
G). Several obscure neutron star subpopulations 
haven been claimed as magnetar candidates, in 
particular Soft Gamma Repeaters (SGRs) and 
Anomalous X-ray Pu ls ars (AXPs); for rev i ews se e 
Woods fc ThompsonI (|2006l ): iMereghettil (|2008l ). 
Most magnetars have been discovered either from 
their persistent X-ray emission properties or when 
they enter into randomly occurring outbursts, dur- 
ing which they emit a multitude of short (^ 100 
milliseconds), soft 7— /hard X— ray bursts. Thus 
far, approximately twenty magnetar sources are 
known and most of them reside in the Galactic 
Plane with a higher concentration close to the cen- 
ter; two are located in the Magellanic Clouds. 

SGRJ0501-K4516was discovered with Swift 
on 2008 August 22, when it emitted a series of 
bright, soft, short bu rsts (JHolland fc Satol l2008t 
Barthelmv et al.l 120081 ) . The first burst also trig- 
gered the Gamma-ray Burst Monitor (GBM) on- 
board the Fermi Gamma-ray Space Telescope 
{FGST; hereafter Fermi). Soon after, our Target 
of Opportunity (ToO) observations with the Rossi 
X-ray Timing Explorer (RXTE) established a pe- 
riod of ~ 5. 76 s in the persisten t X-ray emission of 
the source (JGogiis et al.l 120081 ). Further observa- 
tions with RXTE and the Swift/X-iay Telescope 
(XRT) revealed a spin-down rate of 1.5(5) x 10"^-^ 
s s~^, indic ating a dipole rnagnetic field of 
2.0 X 10^^ G (IWoods et aD l2008l: iRea et al.ll2009l: 
Gogiig et al.l I2OIOI ). Our subsequent Chandra 
ToO observations established an accurate loca- 
tion of the source at R.A.{J2000) = 05h01'"06!76, 
i:»ec(J2000) = -H45°16'33'.^92, w ith an la uncer- 
tainty of O'.'ll (|G6gus et al.ll2010l ). This is the first 
magnetar location at roughly the Galactic anti- 
center direction, placing SGR J0501+4516 most 
likel y at the Perseu s arm of our Galaxy at ^ 2 
kpc (jXuet allbOOel) . 
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The SGR J050H-4516 outburst lasted approxi- 
mately 2 weeks, during which several bright bursts 
were detected with Swift, GBM, RXTE, KONUS- 
Wind, and Suzaku (Enoto et al. 2009; Aptckar 
et al. 2009; Kumar et al. 2010; see also Table 
[1]). After the burst activity ceased, the source 
flux decreased expon entially with an e- fo lding time 
of 27.9 ± 2.5 days (|G5giis et al.ll2010[ ). During 
the entire outburst interval, GBM triggered on 26 
bursts; in addition, an untriggered event search in 
the daily data sets revealed seven more bursts. We 
present here our analyses of the 29 GBM bursts 
from SGR J0501-I-4516 for which we have high- 
resolution data; the properties of the X-ray per- 
sistent so u rce e missi on have been publ ished by 
Rea et al.l (|2009l ) and lGogiis et al.l (|2010t ). In Sec- 



tion |2l we describe the instrument, observations 
and data selection. In Sections |3| and lU we present 
our temporal and spectral results, respectively. 
We discuss the interpretation of our results in Sec- 
tion [51 

2. Instrumentation and Data 

The Fermi/GBM has a wide field of view (8 
sr; un-occulted) and a continuous broad-band en- 
ergy coverage (8 keV - 40 MeV). It consists of 12 
Nal detectors (8 - 1000 keV) arranged in 4 clus- 
ters of three each and 2 BGO detectors (0.2 - 40 
MeV) placed at opp osite sides of the spacecraft 
( Meegan et al.l 120091 ). In trigger mode GBM pro- 



vides three types of science data: CSPEC with 
continuous high spectral resolution (1024 ms and 
128 energy channels), CTIME for continuous high 
time resolution (64 ms and 8 energy channels) , and 
Time- Tagged photon Event (TTE) data (2 /is and 
128 energy channels). For a detailed descrip tion of 



the in strument and data types, see Meega n et al. 
(|2009t) . When a GBM trigger occurs, TTE data 
are provided from ~ 30 s pre-trigger to ~ 300 s 
post-trigger. With its very high temporal resolu- 
tion the TTE datatype is most suitable for the de- 
tailed temporal and spectral analyses of very short 
events like SGR bursts and is, therefore, the only 
datatype used throughout this paper. 

During the burst active period of the source 
(2008 August 22 - 2008 September 3), GBM trig- 
gered on 26 bursts. We further implemented the 
algorithm described in iKaneko et al.l ( 20101 ) to 
search for any untriggered events in the daily data 



sets of 2008 August 21 through 2008 September 
14, and found seven individual short bursts lo- 
cated in the same direction as SGR J0501+4516. 
Among these only three had TTE data and were, 
therefore, included in our analyses. Two of the 29 
bursts were very bright, causing saturation of the 
high speed science data bus; the saturated parts of 
these bursts were excluded from any spectral anal- 
ysis. For all spectra and durations we selected the 
Nal detectors with an angle to the source smaller 
than 50° to avoid attenuation effects. We also ex- 
cluded any detectors blocked by the Fermi/La,rge 
Area Telescope (LAT), or by the spacecraft radi- 
ators or solar panels. The BGO detectors were 
not used as there was no obvious emission in the 
Nals above 200 keV. In columns two through four 
of Table [1] we list for each of the 29 GBM events 
their trigger numbers, times, and the selected Nal 
detectors used for the following analyses. 

3. Temporal Analysis 

Figure [1] exhibits the light curves of four repre- 
sentative SGR J0501-h4516 bursts. Their profiles 
vary from a single short pulse (Figure [1^) to a 
multi-pulse event (Figure [1]:). We note that the 
flat top in the event of Figure [TJl is caused by 
saturation. We have used the GBM TTE data 
to estimate the event durations in both count 
and photon space. Although the former process 
has been used exclusively in th e past to estimate 
SGR and AXP bur st durations (|G6gii!? et al.ll2001 : 
Gavriil et al.l 120041 ) . the latter is used for magne- 
tar bursts for the first time here. We describe the 
methods and the results below; Table [2] contains 
the mean values and widths of the distributions of 
all temporal parameters. 

3.1. Tgo and T^q in count space 

For the Tgdll estima te we used the algorithn i 
originally developed bv iKouveliotou et al.l (|l993l) 
for Gamma Ray Bursts (GRBs) and later adapted 



for S GR bursts (jGogii^ et al.l 120011 : iGavriil et al 



12004^, modified slightly to accommodate the GBM 
TTE data. Each duration was calculated between 
8 ~ 100 keV with 2 ms bins as follows. First we fit 



^Tgo (Tso) is the duration during which the background- 
subtracted cumulative counts increase from 5% (25%) to 
95% (75%) of the total counts. 



the burst background using two time intervals be- 
fore and after each burst ([-2 s, -0.5 s] and [0.5 s, 
2 s]) with a first order polynomial; these intervals 
were kept mostly the same for all bursts unless a 
precursor and/or a tail were detected. Then we 
fit the background-subtracted, cumulative burst 
counts, with a linear plus a step function simulta- 
neously; the linear part was user-selected (before 
and after the burst), while the step part was deter- 
mined from the fit linear trend. The height of the 
step function was then used to represent the net 
total counts (TV) of the burst, subsequently used 
for the Tgo determination. 

Panels a and b in Figure [2] show the distribu- 
tions of Tgo and T50 for all 29 bursts; all individ- 
ual Tgo values can be found in Table [1] (column 4). 
For comparison reasons with other magnetar dura- 
tion distributions (see section 3.4) we fit each dis- 
tribution with log-normal functions and obtained 
(Tgo) = 122.6l^J ms {a = 0.35 ± 0.03, where a is 
the width of the distribution in the log- frame) and 
(T50) = 31.6^2^ ms {a = 0.30±0.03). The average 
values of the raw data weighted by their errors are 
(Tg-o) = 138.3+10°^ ms, and (T^,) = 32Atrs ms. 

The distributions above did not account for the 
uncertainties in the durations. We estimate the 
(max, min) errors (ATgoMAX, ATgoMw) for each 
Tgo (similarly T50) as: 



'Q.95{N+^/N) ^0.05(N~VN) 



TgoMAX = ^o.g5(Ar+VlV) ~ ^1 

TgoMIN = to.95(N-VN) ~ ^Q.05{N+VN) 

ATgoMAX = TgoMAX — Tgo 

ATgoMIN = Tgo — TgoMIN 



(1) 



where i„ is the time when the cumulative light 
curve reaches n counts. We then created proba- 
bihty distribution funct ions (PDFs) for each Tgo 



and Tgo, as describ ed bv lStarling et al.l (|2008[ ) and 



lEvans et al.l (2000). Each PDF was constructed 
using a two-sided normal distribution, where the 
width of each half (a) was set to the uncertainty of 
the duration (ATgoMAX, ATgoMiN, AT50MAX, and 

ATsomin): 



P{x\x,ai,a2) 



V2 



-A 



e^ {x < x) 



y/TT{(Ti +a2) e" {x> x) 

... . A = -{x~xf/2al , , 
withj ; -(2/0 2 (2) 

'^ B = -[X ~ x) /2a2 



Each PDF describes the likelihood to obtain x 
(i.e., Tgo or T50) given its measured value x. Fi- 
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Fig. 1. — Light curves of four bursts from SGR J0501+4516 integrated with 4 ms bin size over 8 — 200 keV. 



nally, we averaged aU sample PDFs to create the 
total PDF of Tgo (T^o) shown in Figure [2] panels 
(a) and (b). Note that the shorter (fainter) events 
have larger errors, resulting in "pulling" the PDF 
towards shorter durations. 
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3.2. Tgg and T^q in photon space 

The photon-based durations, Tqq , are esti- 
mated with an algorithm similar to the one used 
above over each burst cumulative fluence in erg 
cm~^. We used the same time resolution (2 ms) 
and energy range (8 — 100 keV) as in the count du- 
rations. The essential difference here, is that these 
measurements utilize the intrinsic (deconvolved) 
burst spectra instead of the detector recorded 
counts to define the burst intrinsic durations in- 
dependent of different instruments. To perform 
these estimates, we used the GBM public software 
tool R MFIT vSM (for a d escription of this tool 
see also lKaneko et al.ll2006[ ) and the new datatype 
CTTE specially created to facilitate analyses of 
short events. This datatype simply bins the 128 



TTE energy channels into the same 8 bins as the 
CTIME data. The errors in the dura tion estimates 
are tak en from iKoshutJ ()l996al ) and lKoshut et al 
(Il996bl) . 



A detailed description of the photon-based du- 
rations can be found in the Fir st Two Years 
GRB Catalog of the Fermi/GBM ()Paciesas et al. 



in preparationl) . In short, an adequate background 



^http://fermi. gsfc.nasa.gov/ssc/data/analysis/user/ 



interval is selected before and after each burst and 
fit with the lowest acceptable order of a polyno- 
mial to determine the background model parame- 
ters. Next the entire burst interval is fit to deter- 
mine the default set of photon model parameters. 
The model used in these fits is a power law with 
an exponential cutoff (COMPT; described in de- 
tail in Section 4). When all background and source 
model selections are determined for each 2 ms time 
bin, we subtract the background, fit its spectrum 
using the COMPT model, and calculate its pho- 
ton flux. These values are then used as inputs for 
the Tqq {T^q ) estimates, performed with the same 
algorithm described above. 

Figure [3] shows the distributions of Tg^ (T^q ) 
fit with a log- normal function (panels a and b). 
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Fig. 2. — Distributions of T90 (a), T50 (b), rgo (^cj, r5o (^cfj, (59o fej, (^50 (f). The solid lines in panels (a) 
and (b) show the probability distribution functions of T90 and T50, respectively. The dashed curves show 
the best fits with log-normal (c, d) or normal (e, f) distributions. The vertical dot-dashed lines indicate the 
mean values of the fits for each histogram distribution. 



obtaining (T^^) = 124.2^11^ ms (a = 0.38 ±0.06, 
where a is the width of the distribution in the 
log-frame) and (T^o) = 27.6^1;? ms {a = 0.21 ± 
0.03) . The average values of the raw data weighted 
by their errors are (Tj'g 



161-211:6 



and 



50 



49.2 



ms, 
+°-^ ms. The individual T^n values 



-0.8 



■^90 



can be found in Table [T] (column 5) 

The solid curves in Figure [3] exhibit the PDFs 
for TgQ and T^q after taking into account the sym- 
metrized errors in the values. This plot is an 
adaptive kernel density estimation: it was made 
by adding up a set of normalized Gaussian func- 
tions, one at each of the data points and with a la 
width given by the corresponding error estimate. 
The multiple narrow spikes in this figure, mostly 
at the long-duration end of the plot, are largely 
due to the fact that the corresponding errors (in 
contrast to those of shorter and fainter durations) 
are relatively small, yielding tall and narrow Gaus- 
sian components. This is a rather high- variance 
estimate of the distribution function. As an al- 
ternative we show in Figure |4] a completely dif- 
ferent density estimate based on an adaptation of 
the Bayesian block algorithm for histogramming 
of the TgQ and T^q data, taking into account 
measurement errors in the independent variabl e 
(jScarglelllQQSHScargle et aLll2011 in preparation! ). 



This algorithm finds the optimal piece-wise con- 
stant model to represent the data; the optimiza- 
tion corresponds to the maximum likelihood for a 
constant-rate Poisson model for the data in each 
of the bins. The sizes and locations of the bins 
are all determined by this optimization, not pre- 
defined as in ordinary histograms. The solid curve 
is the Bayesian block representation superimposed 
on such an ordinary histogram with bins chosen 
so that the value plotted is more than one only if 
there are duplicate values. This analysis suggests 
that the spiky structure of the PDFs in Figure [3] 
is mostly due to noise fluctuations (the spikes cor- 
responding to a small number of points with small 
formal measurement errors) and does not support 
a quantized or a multimodal duration distribution 
for the bursts from SGR J050H-4516. 

We compare the two duration distributions (in 
count and photon space) in section [ 
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Fig. 3. — Distributions and probability distribu- 
tion functions of Tg'o (a) and TlJ, (b). The his- 
tograms show the raw data and the solid curves 
show the PDFs. The vertical dot-dashed lines in- 
dicate the mean value of the log-normal fits of the 
histograms. 

3.3. T90 (t5o) and (Jgo (i^so) in count space 

The emission time rgg (T50) for each burst was 
determined by adding the time bins (2 ms each) 
of high fluence in decreasing fluence rank un- 
til 90% (50%) of the fluence was reached. The 
emission time interval, thus, characterizes the du- 
ration of high-fluence emission. This parame- 
ter was also first introduced for GRBs as com- 
plementary to their rgn (^50) duration measures 
( Mitrofanov et al.lll999l ). Panels c and d in Figure 
[2] show the distributions of rgo and t^q. These 
are also fit with log-normal functions obtaining 
(t9o) = 70.3+^:5 ms {a = 0.39 ± 0.04) and (rgo) = 




|°gio Tg, 



Fig. 4. — Bayesian block representation of the 
TgQ {lejfj, and T^q {right), taking into account 
measurement errors. The histograms show the 
raw data and the sohd curves show the Bayesian 
blocks. 



20.91:2:3 ms ((7 = 0.30 ±0.05). 

The ratio S90 = rgo/Tgo (Jso = T50/T50) 
is defined as th e duty cycle of each burst by 
Mitrofanov et al.l (|l999l ). This value should not 
be over 1, because the emission time excludes 
low fiuence intervals in the burst, which some- 
times are included in the duration. Panels e 
and f in Figure [5] present the distributions of Sqq 
and Jso, which were fit with normal distributions 
with (^90) = 0.68 ± 0.03 (ct = 0.14 ± 0.03) and 
(^5o) = 0.68 ± 0.02 {a = 0.19 ± 0.02). Both duty 
cycle distributions have the same mean value of 
0.68. 

3.4. Comparisons of durations 

The durations in photon space are estimated 
using the deconvolved source spectra with the 
response of the GBM detectors. In principle, 
since the spectral energies where SGR bursts emit 
most of their photons are relatively narrow (~ 
8 — 100 keV) there should not be a large differ- 
ence between photon and count durations. Fig- 
ure [5] shows Tqq (TJ'q ) versus TqqIT^o)- We notice 
that in general the count space Tqq values tend 
to be larger by a very small amount, mainly due 
to the fact that they take into account number of 
counts irrespective of their energy content. Since 



SGRs have mostly soft spectra, the counts cor- 
responding to the lowest energy photons often do 
not contribute much flux in the durations (i.e., less 
than 5%). The T50 estimates, however, are per- 
fectly aligned along the x = y line, making these 
spectrally-independent measures for durations (as 
also noted for GRBs by Bissaldi et al. 2011). It 
is also important to note here that these photon 
durations validate earlier duration measurements, 
which have been done in count space for all mag- 
netar candidates. 
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Fig. 5. — Toy 'panel: The comparison of the pho- 
ton space rJJ) and the count space Tgo. Bottom 
panel: the similar plot for T50. The solid lines are 
x = y. 

We now compare in Figure [6] the mean values of 
the count space Tgo and rgo of SGR J0501-I-4516, 
with those of four other magnetar candidates. 
SGR J1550 - 4518 was also observed with GBM 

and the temporal parameters are es timated using 

exactly the same pro cedures as here (|van der Horst et al, 
2011 inpreparationi) . SGRs 1806 - 20 and 1900 -h 
14 were estimated us hig RXTE/PCA o bserva- 
tions in 2 - 60keV bv lGdgiis et"aH (l200lh. AXP 
1E2259 -I- 586 durations are from iGavriil et al.l 



(|2004[ ) and are also estimated using RXTE/PCA 
observations in 2 — 60 keV. It is obvious from the 



figure that all durations fall well within the same 
order of magnitude (~ 100 — 150 ms), indicating 
a similar origin for the bursts across the magnetar 
population. 
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Fig. 6. — Mean values of count space Tqq (right 
hatched bars) of four SGRs and one AXP and of 
tqo (left hatched bars) of four SGRs. The photon 
space Tgo of SGR J0501+4516 is indicated with a 
blank bar. Values are marked in ms above each 
bar. 



3.5. 



Relative timing between burst peaks 
and persistent emission pulse phase 



Usi ng the pulse ephemeris reported in lGogiis et al.l 
pOld ). we have aligned each burst peak to the 
pulse phase of the spin period of SGR J050H-4516, 
to search for possible correlations of the burst 
activity with rotational phase. In an effort to 
confirm the GBM barycentric correction for this 
analysis, we compared the barycenter corrected 
light curves of burst bn080823.020 (Table 1) for 
the GBM and RXTE/FCA data. Note that 
the pulse ephemeris is based primarily upon 
RXTE/PCA data. A cross-correlation of the 
GBM and RXTE/FCA burst time histories in- 
dicates no significant shift with an upper limit of 
4 ms. Given the 5.76 s pulse period of this SGR, 
even a 4 ms shift is negligible. As a final check, the 
barycentering software was also tested by epoch 
folding TTE data to obtain a pulse profile of the 
Crab pulsar using the Jodrell Bank ephemeris 
(http://www.jb.man.ac.uk/~pulsar/crab.html). 
We found that the phase of the first peak agreed 
within 200 /is of the RXTE result shown by Rots 



et al. (2004). 

Burst peak phases were computed using the 
event times recorded in the TTE data in the en- 
ergy range 8—60 keV. For each event time his- 
tory, the burst peak time was defined as the av- 
erage event time for the six most closely spaced 
counts in a 4 sec interval surrounding the trig- 
ger time. The phase of pulse maximum was de- 
fined by fitting an inverted parabola to the folded 
RXTE/FCA 2-10 keV light curve. Figure[7]shows 
the distribution of the phase offsets of all burst 
peaks relative to the pulse maximum. The aver- 
age offset (average of the absolute value) is 0.285 
cycles. For a random phase distribution (i.e., null 
hypothesis), one would expect an average offset 
of 0.25 cycles. A Monte Carlo simulation of 10^ 
realizations for 29 draws from a random distribu- 
tion shows that the probability of getting an offset 
value of 0.285 for 29 samples is ~ 26% or roughly 
la. We conclude that there is no correlation of the 
SGR J0501+4516 burst peaks with pulse phase. 
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Fig. 7. — The distribution of the phase offsets of 
all burst peaks relative to the pulsar phase. 



4. Spectral Analysis 

We performed all spectral analysis using RM- 
FIT v3.3 and we generated response files for each 
burst using the GBM response generator gbmrsp 
vl.9; Table [T] lists the detectors used for each 
event. Each event spectrum was restricted within 
8 — 200 keV, as we did not detect any counts above 
200 keV. However, roughly one third of the bursts 
have significant emission in the 150 — 200 keV 



band: the average count rate in the brightest 
detector over the Tgo interval was larger than 
1 count/s. Most of the bursts were faint and they 
only had good statistics for time-integrated spec- 
tral analysis described below. Five events (includ- 
ing the two saturated bursts mentioned in Section 
[2]) were very bright and we were able to perform 
time-resolved analysis as described in Section |4?2] 
Finally, to account for the Iodine K-edge effects 
at 33.2 keV, we excluded from our spectral fits the 
region between 30 — 40 keV. 

4.1. Time-integrated spectra 

We fitted the time-integrated spectra of all 29 
bursts with several models: a single power law 
(SPL), an Optically Thin Thermal Bremsstrahlung 
(OTTB), a single black body (BB), a power law 
with an exponential cut-off (COMPT) , a two 
black body spectrum (BB-I-BB), a single black 
body with a power law (BB-I-PL), a single black 
body with an OTTB (BB-I-OTTB), and finally, 
a single black body with an exponential cut-oflF 
power law (BB+COMPT). The parameters of 
the last two models (BB-HOTTB, BB-hCOMPT) 
could not be constrained by most burst data, while 
a SPL was always a bad fit; all three models were, 
therefore, rejected from further spectral analysis. 
To determine the goodness of fit for the remaining 
models we used the Castor modified Cash-statistic 
(C-stat). This is a modified maximum likelihood 
estimator which asymptotes to x^ , used when 
there are small numbers of counts/bin (Poisson 
regime), which is the case for most of the SGR 
events (especially in the higher energy bins) . 

Using C-stat we were able to further reject sin- 
gle BB, OTTB and BB-hPL models. The first 
two models fit only part of the weaker burst set; 
for these the COMPT model also did not give a 
significantly smaller C-stat. It was not excluded, 
however, because this model, contrary to the first 
two, fit all bursts. Moreover, a COMPT index 
of 1 or —1 reproduces the BB or OTTB spectral 
shape. The BB-fPL model had overall worse C- 



stat values compared to the remaining two models 
(COMPT and BB-^BB). The relative goodness of 
fit among these models is exhibited in Figure [HJ 
which shows spectral fits with COMPT, BB-hBB, 
OTTB, and a single BB of one bright burst from 
SGRJG5G1+4516 (bn080826.136). From the fig- 
ure, we see that the COMPT and BB-hBB models 
can fit the data equally well; the residuals, how- 
ever, of the OTTB and the BB model fits are un- 
acceptably large. 

The COMPT model fits aU 29 bursts weh, with 
the BB+BB model giving equally good fits in 
only 18 events, where we have enough statistics 
to constrain the model parameters. To determine 
whether the COMPT or the BB-^BB model fit the 
data best, we simulated (using RMFIT) a large 
set of bursts with different intensity and spectral 
shape parameters, using COMPT and BB -I- BB 
as input models, and then fit them with both the 
COMPT and BB-^BB models (see van der Horst 
et al. 2011 in preparation, for a detailed descrip- 
tion of the simulations) . We show that the C-stat 
improvement is not significant to conclude that 
BB+BB (with one more parameter) is better than 
COMPT. A set of simulations using Xspec gave 
similar results with RMFIT. Below we discuss our 
COMPT model fits for all 29 bursts. Since the 
BB-f BB fits provid e significant information on the 
source parameters ( Olive et al.l l2004t llsrael et al 



^The analytic expression for this model is: 

/ = Aexp[-S(2 + A)/£;peak](S/Spiv)^, 

where / is the photon number flux in photons s~^ cm~^ 
keV~^, A is the amplitude in the same units as /, -Bpeak is 
the peak energy in keV, A is the photon index, and Epiv = 
20keV is the pivot energy. 



i2008i) , we also describe the results of these fits for 
18 bursts. Section [S] expands on the importance 
of each model. 

11. 1. COMPT model fits 

All SGR J0501-I-4516 bursts were well fit with 
the COMPT model. We list the model parame- 
ters and statistics in Table [1] (columns 6-10). The 
distributions of index and i?poak are displayed in 
Figure [9l The top panel shows the spectral index 
distribution, which is centered around zero and 
is best fit with a normal function with an aver- 
age of -0.32 ± 0.11 {a = 0.9 ± 0.1). This index 
distribution clearly excludes a pure OTTB or BB 
fit for all bursts, as such fits would require in- 
dices of ~ —1, ~ -t-1, respectively. The bottom 
panel of Figure [9] exhibits the Epcaw distribution, 
which also follows a normal function with mean at 
39.8±0.9 keV {a = 9.0 ± 1.0 keV). In both pan- 
els, the hatched areas highlight the distribution of 
18 bursts also fit with the BB+BB model. This 






Fig. 8. — The spectrum of a bright 

SGRJ0501+4516 burst (bn080826.136) fit with 
COMPT, BB+BB, OTTB, and BB models (from 
top to bottom). 



subsample was fit with a normal distribution with 
(i;peak) = 36.5 ± 1.6 keV, ct = 6.1 ± 1.4 keV and 
(index) = -0.63 ± 0.04, ct = 0.62 ± 0.04. Most 
bursts here have negative photon index and lower 
Epcak, indicating a softer spectrum. 



ITi 

in 






Z [ 





- 




./■-, 


\ 




- 


/ 






^ 


\ 
\ 

\ 
\ 


\ 

\ 



4 
2 

-3 -2 - 



1 

index A 









1z 


; 


10 


/ 


/ 


\ 


- 


8 


^ 


\ 

\ 


\ 


6 


- 




;^^ 


\ 


\ 


4 


- 


9 


\ 


- 


o 




W//A 


\ 


- 


/ 


^%^i^%^^^^^i%^^^<i^%^^^^^^^^^ 


s 


- 


n 


^^ .■' 


;^^^^^^^^^ 


"-> 


- 



10 20 30 40 50 60 70 
Epeok (keV) 

Fig. 9. — Distributions of index (top panel) and 
Epeak (bottom panel) of the COMPT model fits 
for 29 bursts from SGR J0501+4516. The right 
slashed bars represent the subset of the 18 bursts 
that can be fit with the BB+BB model as well. 
The dashed and dotted lines are the best fits with 
normal distributions of all bursts and of the subset 
of 18, respectively. 

Using the COMPT spectral fits we estimated 
the event fiuences (8 — 200 keV, also listed in Ta- 
ble [Ij and plotted in Figure [10] their correlation 
with spectral indices (top panel) and with ii^pcak 
values (bottom panel) . We note that the brightest 
events have a constant index value of ^ 0, which 
progresses to lower (softer) values for weaker ones 
and is widely scattered with larger errors for the 
faintest events. Similarly, a simple trend can not 



10 



describe the relation between the hardness (Epcak) 
and fluence. For bursts with high energy fluence 
(e.g., > 5 X 10"^ erg cm"^), Speak values are con- 
stant and cluster around 35 keV, while the values 
of weaker bursts cluster around 45 keV with a large 
scatter range between 30 and 60 keV. We note here 
the large span of fluence in our data (3 x 10~^ to 
3 X 10~^ erg cm~^), which is one order of magni- 



tude broader than earlier results (IFenimore et al 



19941 : iGoeus et al.ll200ll: [Cavriil et alJl2004l) . We 



expand on these trends in Section [5j 



^ 






2F 




10"^ 10"^ 10"^ 10" 



100 



> 

0) 



energy fluence (erg/cm^ 






10" 



0" 



0" 



0" 



energy fluence (erg/cm^ 



Fig. 10.— Scatter plots of the COMPT index {top 
panel) and -Epeak {bottom panel) versus the event 
energy fluence (8 — 200 keV) for 29 bursts from 
SGRJ0501+4516. 

We now proceed to estimate the cumula- 
tive energy fluence during the active period of 
SGR J050H-4516, shown at the top panel of Fig- 
ure [11] The plot levels off at 3.75 x 10"^ erg 
cm~^, which should be considered a lower limit 
(as we have not taken into account event sat- 
uration, untriggered events without TTE data. 



and additional bursts seen, e.g., with Swift). As- 
suming that the source is 2 kpc away from the 
Earth, then this fluence would correspond to a 
total energy of at least 1.8 xlO^" erg emitted 
from the magnetar in bursts during this 13-day 
active period. The bottom panel in Figure [11] 
presents the differential distribution of energy flu- 
ence {\ogN—\ogS). The best fit with a power law 
function is overplotted in the figure. The index 
of the power law is —0.48 ± 0.02, which corre- 
sponds to dN/dF oc F-i-48_ This slope is similar 
to the one es timated with the RXTE/VGPl for 
SGR 1806-20 (|Gogiis et al.ll200fl[ ). but differs from 
all other SGR and AXP sl ope estimates, which 
are a ll very close to —1.7 ([Woods fc Thompson 
20061 ). It also differs from the slope estimate of 
the SGR 1806 - 20 events detected with the Inter- 
national Cometary Explorer (ICE) and the Burst 
And Transient Source Experiment (BATS E) data, 
which is -1.7 as well ([Gogiis et al.|[2000l) . How- 
ever, the GBM fluence range covers the higher end 
of fluence values, and the analyzed burst sample 
is most likely incomplete at the lower fluence end 
(which contributes most of the events in all other 
sources) . 

4.1.2. BB+BB model fits 

Eighteen bursts were spectrally fit equally well 
with the BB+BB model. We chose to also present 
these results here, although this model is more 
complicated than the COMPT model (with one 
more parameter), because it has been success- 
fully used in the past for the s pectral analysis of 
SGR bursts ([Oliye et aL|[2004 fperoci et al.l 12004 : 



[Israel et al.l l2008^: we will compare our spectral 
results in Section [5] Assuming that the two BB 
components arise from two hot spots on the sur- 
face or photosphere of the source, each emission 
area, i?^, can be calculated from the temperature, 
T, of the corresponding black body spectrum as: 



i?2 = FD^/aT^ 



(3) 



where F is the average energy fiux per event (total 
burst fluence divided by the spectral integration 
time) , D is the distance to the magnetar (assumed 
to be 2 kpc), and a is the Stefan Boltzmann con- 
stant. 

We plot the emission area of both BB compo- 
nents as a function of fc T in Figure [T^ As shown 
earlier by [Israel et al.[ ([2008[) for SGR1900 + 14, 
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Fig. 11. — to'p 'panel: The evolution of the inte- 
grated energy fluence of SGR J0501-I-4516, calcu- 
lated with a COMPT spectral model. The cumu- 
lative plot starts at the first trigger on 2008 August 
22. bottom panel: The differential distribution of 
energy fluence estimated with a COMPT model. 
The dashed line is the best fit with a single power 
law of index -0.48 ± 0.02. 



we also note a clear separation of the two temper- 
atures and emission areas of the cooler and hotter 
black bodies. The cooler BB has a larger emission 
area, and the temperature spreads between 3 — 7 
keV. The emission area of the hotter BB with tem- 
peratures in the 10 — 20 keV range is much smaller. 
The emission areas of both BB components have 
similar evolution through time during the active 
period of SGR J05014-4516 (Figure [l3l). The rel- 
ative increase in emission area, however, is higher 
(by a factor of ten) for the hotter BB component. 
Finally, we notice in Figure [14] that the total flu- 
ence of the burst is divided equally between the 
two black body components (we find a correlation 
coefficient of 0.94 corresponding to a chance prob- 



ability of 5.21 X 10 ). A single power-law with 
an index of 1.00 ± 0.05 fits the data well. 
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Fig. 12. — Emission area as a function of black 
body temperature for time-integrated spectra. 
The dots mark the black body component with the 
higher temperature, while the circles represent the 
lower temperature black body. The upward arrows 
indicate the saturated bursts. 
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indicate the saturated bursts. 



4.2. Time-resolved spectra 

We performed time-resolved spectral analysis 
with 8 ms temporal resolution for the five brightest 
bursts of SGRJ050H-4516. We binned the data 
requiring a significance of at least 3a above back- 
ground for each bin. Each bin was then fit with the 
COMPT model, which was the best model for the 
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Fig. 14. — The correlation between time- 
integrated fluence of cooler BB and hotter BB. 



time-integrated spectra. Figure [15] displays the 
light curves of these five events overplotted with 
their -Epeak values (left column) and the correla- 
tion between -Epeak and energy flux (right column) . 
We see that the Ep^ak follows the light curve for 
the brightest part of the burst in four out of five 
cases (except for the fourth panel from the top), 
and it rises surprisingly at the beginning and the 
tail end of each event. 

To investigate this trend, we plotted in Figure 
[in] the combined -Epoak values from all five bursts 
versus their fluxes (top panel) . We clearly see that 
E^pcak rises at both high and low flux values with a 
minimum determined with a broken power-law fit 
at 8.7±0.9x lO^^ergcm^^ s^^. The same trend is 
seen in the lower panel between i?poak and fluence 
with a minimum at 3.4±0.2 x 10~^erg cm~^. The 
indices below and above this critical flux (fiuence) 
values are -0.28±0.03 and 0.14±0.01 (-0.14±0.01 
and 0.28 ± 0.02), respectively. An F-test compar- 
ison of a broken power law and a single power 
law fit shows that the former provides a better 
description of the correlation, with a small prob- 
ability of chance coincidence 1.2 x 10~^ for the 
flux and 9.8 x 10^^ for the fluence. We also per- 
formed a Spearman rank test for both branches of 
each plot and found a low flux (fluence) correlation 
coefficient of —0.66 (—0.59) with a probability of 
chance coincidence of 1.7x10"^ (2.7x10""). The 
correlation coefficient and probability for the high 



flux (fluence) parts are 0.51 (0.78) and 8.7 x lO"*' 
(3.3 X 10~^), respectively. The latter probabili- 
ties are less significant than the ones for the low 
flux/fluence parts, because the Spearman rank 
tests involve a smaller number of data points. For 
the full flux (fluence) data sets, the correlations are 
not very signiflcant, with a coefficient of —0.23 (- 
0.28) and probability of 6.4x10-3 (1.3x10-^), giv- 
ing further support to changing trends (minima) 
in the correlations between -Epoak and flux/fluence. 
Figure [17] displays the distribution of the 
COMPT index and Epcak of all time-resolved spec- 
tra. Both distributions were fit with normal func- 
tions with (Epeak) = 33.7 ±0.5 keV {a = 8.5 ±0.5 
keV) and (index) = -0.11±0.09 {a = 0.47±0.09). 
These values are very similar to those obtained 
with the time-integrated fits. 

Finally, we fit the time-resolved spectra with a 
BB+BB model. Retaining only the fits which con- 
strained the BB+BB parameters, we calculated 
the emission areas for each BB component. Simi- 
lar to the integrated spectra. Figure [18] shows that 
the emission areas follow different behaviors with 
temperature for the hard and soft BB. Interest- 
ingly, a co mparison with sim ilar results shown in 
Figure 5 of (" Israel et al.ll2008i) for SGR1900+14 re- 
veals that although the kT values are very close for 
the two BBs, the emission areas differ by at least 
one (and maybe two) orders of magnitude. These 
results indicate that either the two sources have 
different burst energetics or that (very unlikely) 
the distance determination of SGR1900 + 14 is off 
by a large factor. 

Using the 2 kpc distance for the source, we 
estimate the isotropic luminosity of the two BB 
components shown in Figure 1191 The Spearman 
rank correlation coefficient is 0.83 corresponding 
to a chance coincidence probability of2.78xl0-^^, 
indicating that these components are well corre- 
lated, as also shown with a powe r-law fit with in- 
dex 1.1 ± 0.1. llsrael et al.l ( 20081 ) have performed 
a similar analysis for the data of SGR1900 + 14 
and found a spectral index of 0.70 ± 0.03. 

5. Discussion 

5.1. COMPT versus BB+BB 

Our spectroscopic fitting clearly indicates that 
the COMPT and BB+BB models yield superior 
fits to the other possibilities. The COMPT model, 
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Fig. 15. — Left column: Light curves of the five brightest events from SGR J0501+4516, overplotted with 
their ii^peak values. The dark areas mark the saturated parts in two of the bursts. Right column: Correlation 

between £^peak and energy flux in the 8 — 200 keV band. 
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Fig. 16. — (top panel): Correlation between the 
-E'peak and energy flux in the time-resolved spectra 
of the five brightest bursts of SGR J0501+4516. 
[bottom panel): Correlation between the -Epeak 
and energy fluence for the same bursts. The solid 
lines exhibit the best fit to the data with broken 
power-law functions. 
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Fig. 17. — Distribution of index (top panel) and 
-E'peak {bottom panel) of the COMPT model for all 
time-resolved spectra. The dashed lines are the 
best fits with a normal distribution. 



with its power-law shape curtailed by an expo- 
nential turnover, is intended to mock-up the clas- 
sic unsaturated Comptonization spectrum real- 
ized in models of accretion disks such as in Cyg 
X-1 or in active galactic nuclei (see Chapter 7 
of Rybicki & Lightman 1979 for a summary of 
its development as a solution of the Kompaneets 
equation). These models use hot, thermal elec- 
trons in a corona to repeatedly scatter low-energy 
photons, heating them gradually up to an en- 
ergy E ^ kTe consistent with the electron tem- 
perature Te, at which point a quasi-exponential 
spectral turnover emerges as further heating be- 
comes impossible. The power law marks the scale- 
independence of the Compton upscattering, and 
its slope depends only on the mean energy gain per 
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Fig. 18. — Emission area as a function of black 
body temperature for time-resolved spectra. The 
dots mark the black body component with the 
higher temperature, while the circles represent the 
lower temperature black body. 

collision ((A£') = 4fcTe for non-relativistic elec- 
trons) and the probability of loss of photons from 
the scattering zone, being expressed as 



dN 
dE 



ex E'' 



-A = — 



9 4 
4 y 



(4) 



with the Compton y-parameter in the domain 
y < 1. This parameter is the product of the 
average fractional energy change per scatter- 
ing and the mean number of scatterings, and is 
4fcTe/(TOeC^) max{T, r^} for a non-relativistic sit- 
uation for repeated Compton upscattering. Here 
T is the scattering optical depth. 

In the context of magnetars, such Comptoniza- 
tion can also ensue, but the strong magnetic field 
now plays an important role. Ephemeral coronae 
of hot electrons could be expected in a dynamic 
inner magnetosphere. For example, this could be 
due to intense dissipation of magnetic energy in 
the closed field line region via field line twisting, 
i.e., transient departures from poloidal geometry, 
as in the considerations of Thomson, Lyutikov 
& Kulkarni (2002), Thompson & Beloborodov 
(2005), Beloborodov & Thompson (2007), and No- 
bili, Turolla & Zane (2008). Such a pumping 
of energy into electrons in low altitude regions 
would then be subject to irradiation by the in- 
tense bath of surface X-ray emission. The elec- 
tron coronae would mimic those of their black 
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Fig. 19. — The correlation between the time- 
resolved luminosities (in 10'*° erg s^^) of the soft 
and hard BE. 



hole counterparts discussed above, and serve as 
a Comptonization target for the X-rays. Tempo- 
rally, the coronae could be quite variable, resulting 
in varying or chaotic time profiles such as are ob- 
served. If the deposition of energy in hot electrons 
is persistent over many light-crossing timescales, 
so also can the upscattered hard X-ray emission 
be. In this magnetic case, the spectrum would 
again naturally be a power law truncated by an 
exponential tail, whose energy is pinned by T^. 
The slope does differ somewhat in value from the 
non-magnetic case because the presence of the in- 
tense field anisotropizes the environment. This al- 
ters the average fractional energy change per scat- 
tering from the isotropic 4fcTe form in Rybicki & 
Lightman (1979). We note that this picture is sim- 
ilar to that of Lyutikov & Gavriil (2006; see also 
Rea et al. 2008), who modeled the steep X-ray 
tails below 10 keV in quiescent magnetar emission 
using a resonant cyclotron scattering picture. 

Since the electrons will move along the field 
lines in the zeroth Landau level, the exact kine- 
matics that impact the determination of (AE) de- 
pend on the colatitude and altitude of the col- 
lisions. Therefore, the effective magnetic Comp- 
ton y -parameter, ys, that could be substituted in 
Eq. (j4|), would take a somewhat different value 
from its non-magnetic cousin, but one anticipates 
that the range of spectral realizations would be 
similar to the B — Q case. Since head-on colli- 
sions generally yield greater heating of photons. 
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and these are more likely at lower altitudes, it is 
expected that interaction zones close to the stellar 
surface would spawn larger ys and therefore flat- 
ter photon spectra. If there were a coronal radius 
expansion, this would then translate to an evolu- 
tionary steepening of the spectra with time. Note 
that the photon retention probabilities in the coro- 
nae will need to be higher than in the Lyutikov & 
Gavriil (2006) scenario to generate the requisite 
moderate optical depths r to match the flat out- 
burst spectra discussed in this paper. 

The scatterings will take place below the cy- 
clotron resonance unless their altitude is above 
around 10 stellar radii (e.g. Lyutikov & Gavriil 
2006). Below resonance, the cross section is far in- 
ferior to the Thomson value (e.g. see Herold 1979). 
Then, the coronal electron density must accord- 
ingly be much higher than when i? = in order 
to establish a sizable optical depth. This provides 
a possible distinction between the two classes of 
magnetar emission: the steady < lOkeV signal 
may originate at higher altitudes where the reso- 
nance is accessed and the local electron density is 
low, whereas the bursts reported here may be trig- 
gered closer to the surface in higher density zones 
that initially precipitate scattering below the cy- 
clotron resonance. As photons get energized above 
20 keV, the altitude where the resonance is ac- 
cessed is lowered, increasing the cross section for 
scattering and therefore ys- Note also, that to a 
considerable extent, polarization mode-switching 
(e.g. Miller 1995) between the two photon po- 
larizations can help increase the opacity, adding 
a further nuance to consider when modeling the 
scattering environment. 

Observationally, it is difficult to discern unam- 
biguously the presence of the intense field using 
a Comptonization model scenario in the energy 
range of data presented in this paper: details of 
both the field strength and the geometry are sub- 
sumed in a single parameter ys- Since the emis- 
sion in this scenario should be strongly polarized, 
and the degree of polarization should depend on 
the interaction geometry, a hard X-ray polarimeter 
would provide insightful probes into the presence 
of a strongly anisotropizing super-critical field. 

The dual blackbody can also be envisaged as 
a viable alternative from a theoretical standpoint. 
The moderate Thomson depths required to gener- 
ate fiat Comptonization spectra could easily be 



higher, thereby pushing the electron-photon in- 
teractions more towards an equilibration. The 
saturation temperature is then controlled by the 
largely uncertain total energy dissipated in the in- 
ner magnetosphere per hot electron present. Since 
any equilibration will be non-uniform over a coro- 
nal volume, there should be a modest tempera- 
ture gradient throughout, smearing out the con- 
tinuum. The range of temperatures will not be 
great because the system is not gravitationally hy- 
drostatic in character. Moreover, radiative trans- 
fer effects impact the spectral shape and further 
modify it. Accordingly, pure blackbody shapes are 
not expected. It is quite conceivable that a two- 
component blackbody fit may well approximate 
the emergent continuum that is a superposition 
of distorted blackbodies spanning a small range of 
temperatures. Most probably, due to general ther- 
modynamic considerations, the base of the coro- 
nae (e.g. centered near the source of magnetic 
dissipation) should be hotter than the outer layers 
(see also Lyubarsky 2002 and Thompson & Dun- 
can 1995). Interestingly, the fits here generate a 
smaller volume associated with the hotter black- 
body contribution, consistent with the expectation 
for coronal structure. Yet, somehow, we must ob- 
tain a view of the hotter zone, therefore indicat- 
ing a strongly aspherical coronal geometry. Mod- 
eling this semi-equilibration is a challenging task 
for theorists considering the influences of the field, 
the twisted magnetospheric geometry, and the in- 
herent anisotropy and polarization-dependence of 
the scattering process. At present, it is not pos- 
sible to distinguish between this thermal scenario, 
and a Comptonization one. 

5.2. Conclusions 

Since the COMPT model fits all events, we 
used it to derive fluences (time-integrated and 
time-resolved), and determined - using the time- 
resolved values - that the hardness-fluence corre- 
lation can be described by a broken power law, 
with a minimum at a fluence of 1.7 x 10^^ erg 
cm~^. We have used £^poak values to charac- 
terize the hardness of the events in our sam- 
ple. Earlier studies of this correlation used hard- 
ness ratios to bypass proper spectral fits due 
to low qualit y or insufficient data . Moreover, 



these studies (jFenimore et al.l 119941 : iGogiis et al. 



20011 ) had a very narrow overlap in fluence space. 
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Fenimore et aD (Il994l) used 95 SGR1806 - 20 
events detected with the Interplanetary Sun Earth 
Explorer-3 (ISEE-3) with fluences ranging be- 
tween 1.25 X 10^^ — 8 X 10^^ erg cni^^, and 
found a shghtly positive correlation between the 
two quantities (hardness increasing with fluence). 
On the other hand, iGogiis et alJ (|200l[ ) analyzed 
159 and 385 events from SGRs 1806 - 20 and 
1900 + 14, respectively, observed with RXTE, with 
fluences ranging between 1.0 x 10^^ — 2.0 x 10^^ 
erg cm~^, and found the opposite trend (an anti- 
correlation between hardness ratios and fluences) . 
Since the GBM data sample covers both fluence 
ranges (2.0 x 10~* — 2.0 x 10~^ erg cm~^) we were 
able to establish that both trends are indeed cor- 
rect, and to define the turning point in the iJpoak 
- fluence diagram. 

Only two magnetar candidates were observed 
with GBM to emit a multitude of bursts, thus 
allowing us to construct their E-p^ak - flux dia- 
grams. We find for the time-resolved data of 
SGR J0501-I-4516, that E'pcak reaches a minimum 
of -- 30keV at a fiux value of 8.7 x lO"*" erg cm"^ 
s"^ The second source (SGRJ1550- 5418; van 
der Horst et al. 2011 in preparation) has fiux 
values of 4.4 x 10~^ erg cm~^ s~^ at a similar 
-Epcak minimum. In addition we used the hardness 
ratio - count rate relationsh ip of SGR 1806 — 20 



found in the Integral data ( Gotz et al.l l2004r ) to 



derive an approximate flux value of 1.7 x 10~^ erg 

cm^^ s^^ (using their conversion from count rate 

to flux) at the hardness ratio minimum (Figure 

I — — 1 1 — — k 

3 in lGotz et al.M2004[ ). We then converted these 

values into isotropic source luminosity, iiso (see 
Table [3]) ; we note that (with the caveat of small 
number statistics, and uncertainties in the dis- 
tance measurements) these values are comparable 
[(0.4-1.5) X 10'*° erg s"!], although the B-fields 
and fluxes at E'peak minima vary by approximately 
a factor of ten between the two GBM sources and 
SGR 1806 - 20. Whether these differences reflect 
intrinsic source properties or instrumental effects 
is not yet clear. The ii^pcak trend is, however, 
clearly established in at least two different instru- 
ment data sets. The physical interpretation of this 
trend is beyond the scope of this paper. 

Part of our sample (18 bursts) was fit with 
the BB-I-BB spectral function applied to time- 
integrated intervals; the remaining 11 events could 
not be fit due to poor statistics. In addition for five 



bright events we performed time-resolved spectral 
analysis. We found that the temperatures, emis- 
sion areas and fluences (luminosities) of the two 
thermal components exhibit very similar proper- 
ties in both cases (time-integrated and -resolved). 
Our results are consistent with those presented for 
int ermediate and sho rt bu rsts from SGR 1900 -|- 14 
bv lOliveet all (|2004l ) and llsrael et~ai] (|2008[ ): the 
SGR J0501-I-4516 burst emission areas and tem- 
peratures fall into the same region with those of 
the shor t bursts of SGR 19 00 -I- 14 shown in Fig- 



ure 5 of llsrael et al.l (|2008[ ). We also see a sim- 
ilar trend in the correlation between the lumi- 
nosities of t he tw o components as described by 



Israel et al.l (|2008l ). namely that both luminosities 
increase in tandem. This behavior indicates that 
the hot and cool BB components may come from 
tw o separate emission regions, as also pointed out 
by llsrael et al.l ([20081) : a smaller but hotter one 
from the surface of the magnetar, and a larger, 
cooler one from the star's magnetosphere (but see 
also the discussion in section 5.1). We note here 
that the time-integrated hot BB emission area 
of 13 SGR J0501-K4516 bursts (i.e., not including 
the 5 brightest events; Figure [T2|) is similar to 
the emission area (~ 0.05 km^) of the BB com- 
ponent found in the persistent emission of SGR 
J1550— 5418 during one active bursting episode in 
January 2009, believed to originat e from a hot 



spot on the neutron star surface (jKaneko et al 
I2OIOI) . 

Finally, we performed a detailed temporal anal- 
ysis of all 29 bursts and estimated their durations 
(T90/T50) for the first time in count and in pho- 
ton space. Both estimates agree within statis- 
tics and thus validated all earlier (count space) 
duration estimates of magnetar candidate bursts. 
The durations (and for four SGRs also the emis- 
sion times and duty cycles) of five magnetar can- 
didates follow a log-normal or normal distribu- 
tion. We find that SGR J0501-I-4516 events are 
very similar in average duration with four more 
magnetar candidates (three SGRs and one AXP). 
However, the Tgo distribution of the bursts from 
AXP 1E2259 + 58 6 has a factor of tw o larger dis- 
persion {a ^ 0.73. iGavriil et al.ll2004l ) than those 
of at least three other SGRs [a ~ 0.34, 0.35, 
Gogiis et al.l I2OOII and ~ 0.35, present work). 



SGR J0501-I-4516 bursts have an average duty cy- 
cle (^90 = 0.68) larger than other SGRs (0.45, 
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0.46 for SGRs 190 + 14, 1806 - 20, respectively; 
Gogtis et al.ll20Ql ) . The differences in the intrinsic 
properties of the sources may be due to differences 
in the size of the active region responsible for the 
burst emission. We discuss below the burst ener- 
getics and its evolution. Similarly to other mag- 
netars, we also do not find a correlation between 
the pulse phase of the source and the burst peak 
times. 

The overall behavior of SGR J0501+4516 during 
its 13-day active period (2008 August 22 to 
September 3) is also very interesting. Although 
about half of the bursts (16/29) occurred on one 
day (2008 August 23), most of the burst energy 
was emitted from the source during August 24 — 26 
(see also Figure [TT|) . The average fluence of the 
beginning and the later part of the active period 
was constant and at a lower level. During these 
three days, GBM detected seven bursts, five of 
which are also distinguished in that the emission 
areas from both BB components are the largest. 
Further, one of these five events (bn080824.054) 
shows a double-peaked structure reminiscent of 
some bright thermonuclear Type I X-ray bursts 
from accreting neutron stars, which has been in- 
terpreted as due to p hotospheric radius expansion 
(PRE). Watts et all ( ^010) studied the effects of 
PRE in high magnetic fields using this event and 
find that the predicted flux from PRE theory is 
consistent with the one observed, opening the way 
to determining fundamental parameters of neu- 
tron stars, such as their equation of state. 
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Table 1 
Summary of time-integrated spectral analysis for the 29 GBM bursts from 

SGRJ0501+4516. 



No. 



Trigger catalog 

# 



Trigger time 
(UT) 



Detectors 



T90 
(ms) 



rppn 
-^90 

(ms) 



Spectral 
index'^ 



(kcV) 



C-stat/doF 



Flu 



Peak flu 



1 


bn080822.529» 


12:41:56.914 


8, 


7, 4 


86t- 


80 ± 16 


1.06 ±0.70 


40.83 ± 2.76 


209.37/177 


7.05 ± 0.62 


1.59 ±0.32 


2 


bn080822.647t 


15:36:35.200 


9, 


10 


216tfo 


226 ± 24 


-1.32 ±0.34 


39.83 ± 5.98 


147.70/116 


19.3 ± 1.42 


7.04 ±0.67 


3 


bn080822.981 


23:32:57.746 


2 




30+_ll 


30 ± 15 


1.48 ± 1.42 


44.86 ± 5.02 


47.19/57 


4.41 ±0.67 


2.37 ±0.52 


4 


bn080823.020''''' 


00:28:09.904 


3, 


4 


661- 


48 ±7 


-1.27 ±0.20 


36.37 ±3.28 


149.25/117 


25.02 ± 1.12 


6.61 ±0.59 


5 


bn080823.091'' 


02:11:36.630 


IC 


), 11 


676tlt 


554 ± 40 


-1.17±0.17 


42.09 ±2.63 


154.94/117 


82.84 ±3.04 


6.49 ±0.74 


6 


bn080823.174 


04:10:19.280 


0, 


1 


M7tll 


330 ± 51 


-0.51 ±0.44 


57.53 ± 7.06 


130.57/119 


14.13 ± 1.29 


1.63 ±0.37 


7 


bn080823.248 


05:56:31.529 


2 




„-„ + 131 


276 ± 34 


1.03 ±0.55 


51.94 ±3.41 


65.11/57 


22.18 ± 1.80 


3.42 ±0.64 


8 


bn080823.293'' 


07:01:09.967 


3, 


0, 1, 5 


174+70 
-^'^-20 


164 ±7 


0.52 ±0.29 


48.13 ±1.85 


272.41/239 


20.10 ±0.90 


2.89 ±0.34 


9 


bn080823.293+ 


07:04:22.610 


3, 


0, 1, 5 


00 + 24 


30 ± 11 


-1.63 ±0.26 


26.68 ±8.25 


276.95/240 


9.54 ±0.56 


5.59 ±0.36 


10 


bn080823.319 


07:39:32.257 


9, 


10 


1421™ 


122 ± 25 


-0.98 ±0.30 


36.96 ±3.22 


14.17/115 


19.42 ± 1.16 


4.03 ±0.49 


11 


bn080823.330 


07:55:45.690 


4, 


3, 8, 7 


1921™ 


162 ± 13 


-0.79 ±0.14 


30.10 ± 1.10 


282.62/238 


67.05 ± 1.55 


15.24 ±0.71 


12 


bn080823.354'' 


08:30:01.633 


11 


-, 8 


96tlf 


94 ± 114 


0.09 ±0.80 


28.10 ±2.89 


141.75/119 


8.62 ±0.83 


2.97 ±0.37 


13 


bn080823.429 


10:18:13.891 


0, 


1, 3, 5 


94t- 


82 ± 13 


-0.85 ±0.24 


55.32 ±4.73 


262.81/238 


14.24 ±0.76 


5.02 ±0.40 


14 


bn080823.478°''^ 


11:27:32.306 


8, 


4 


264if, 


246 ±6 


-0.12 ±0.10 


34.50 ± 0.48 


111.42/118 


512.6 ±6.48 


69.62 ±2.80 


15 


bn080823.623° 


14:56:23.563 


IC 


), 11 


2201^2 


204 ±21 


0.88 ±0.51 


48.70 ±2.89 


128.45/118 


21.12 ± 1.50 


3.02 ±0.54 


16 


bn080823.714 


17:08:49.038 


9, 


10 


406tf« 


398 ± 11 


1.50 ±0.41 


46.68 ±1.69 


133.23/116 


33.04 ± 1.59 


3.11 ±0.47 


17 


bn080823.847 


20:19:30.659 


9, 


10 


264lf3„ 


124 ± 11 


-1.12 ±0.14 


36.57 ± 1.77 


106.81/116 


78.61 ±2.26 


19.95 ± 1.04 


18 


bn080823.847+ 


20:23:42.822 


9, 


10 


,„„ + 224, 
108_52 


110 ± 108 


-1.01 ±0.23 


29.95 ±2.22 


120.36/116 


33.09 ± 1.31 


10.29 ±0.71 


19 


bn080823.986'' 


23:39:24.472 


9, 


11, 7, 6 


60tfe 


30 ± 10 


-0.37 ±0.60 


46.62 ± 5.84 


251.91/241 


4.37 ±0.49 


1.43 ±0.28 


20 


bn080824.054* 


01:17:55.394 


2 


5 


260t« 


250 ±3 


-0.25 ±0.04 


36.01 ±0.24 


216.46/119 


1537 ±9.88 


185.90 ±6.49 


21 


bn080824.346 


08:18:24.418 


3, 


4 


34+«« 


28 ± 11 


-0.22 ±0.48 


57.33 ± 7.0 


107.47/118 


5.70 ±0.61 


3.29 ±0.50 


22 


bn080824.828 


19:52:51.264 


2, 


10 


82t- 


62 ± 16 


0.5 ±0.84 


43.22 ±3.98 


127.93/116 


6.39 ±0.72 


1.64 ±0.36 


23 


bn080825.200'^'* 


04:48:27.405 


4 




not- 


102 ±8 


-0.36 ±0.12 


30.32 ±0.7 


73.1/58 


213.6 ±4.10 


103.07 ±3.70 


24 


bn080825.401 


09:37:42.158 


4, 


3, 8 


128+1^ 


114 ±4 


-0.03 ±0.13 


37.76 ±0.66 


211.69/176 


104.8 ±2.01 


36.99 ± 1.28 


25 


bn080826.136'^''' 


03:16:14.773 


8 




160+^^ 


146 ±7 


-0.09 ±0.1 


36.51 ±0.56 


74.49/58 


507.3 ± 7.78 


161.30 ±5.85 


26 


bn080826.236 


05:40:19.425 


9, 


10 


88t- 


100 ± 16 


0.03 ±0.36 


51.88 ±3.41 


90.78/116 


17.08 ± 1.06 


4.74 ±0.57 


27 


bn080828.875 


20:59:39.966 


1, 


0, 5, 3 


72+_tt 


44 ±22 


-0.85 ± 0.47 


43.41 ±5.69 


245.95/239 


5.28 ±0.51 


1.48 ±0.24 


28 


bn080903.421 


10:06:35.329 


4, 


5 


50+11 


60 ±83 


-1.07 ±0.48 


47.35 ± 7.63 


134.69/119 


10.96 ± 1.09 


4.06 ±0.67 


29 


bn080903.787 


18:53:48.775 


2, 


10 


100+',^ 


80 ±6 


-0.95 ± 0.39 


33.68 ±3.75 


128.86/117 


13.88 ± 1.01 


6.68 ±0.65 



Note. — * saturation burst; ' untriggcrcd burst 
" Simultaneously detected with Swift/HAT 
' Simultaneously detected with RXTE /PCA 
''Simultaneously detected with KONUS /Wind 

Simultaneously detected with Suzaku 
"Calculated with the COMPT model 
•'in 10~*erg cm"^ between 8 - 200 keV 
"In 10~"'erg cm"^ s"^ for 8 ms between 8 - 200 keV 



Table 2 
Parameters of duration distributions and weighted mean durations for 29 burst from 

SGRJ0501+4516. 



T90 



T50 



Mo 



Weighted inean'^ 



122.6t.l'l 
0.35 ±0.03 
138.3t-- 



31.6'' 



0.30 ±0.03 

,+0.9 



32. 4^^ 



0.8 



1 24 Q"*"^^'"^ 
-^^^•^-15.2 

0.38 ±0.06 
161.2 ± 1.6 



27.6t\■^^ 
0.21 ±0.03 
49.2 ±0.8 



70.3i^;^, 
0.39 ±0.04 



-,+2.5 

0.30 ±0.05 



20.9^ 



0.68 ±0.03 0.68 ±0.02 
0.14 ±0.03 0.19 ±0.02 



Note. — '^ ms in columns 1 — 6, dimensionless in eolumns 7 &: 8; in the log-frame except for (5go and d'so; ^ in ms. 



Table 3 
iiso at the hardness turnover for three magnetar candidates. 



Distance 
(kpc) 


Flux 
(IQ-s erg cm-2 g-i 


2a 

3.9" 

8.7'= 


8.7 
4.4 
1.7 



-^iso -^surface 

(10*0 erg s-i) (10" G) 



SGRJ0501+4516 

SGRJ1550-5418 

SGR1806-20 



0.41 


2.0^ 


0.82 


2.2^^ 


1.5 


20.6f 



Note.— °IXu et al.l (l2006h: HGogiis et all (1201(1): '^ iTiengo et all (l2010h : 
Camilo et all (|2007l ): ^Bibbv et al.l (|2008l ): 

■f The mean surface d ipole field ^surface = 3.2 x lO^^vPP G, P and P frc 



Mereghetti et al.l (|2005l ). 
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